QCD Phase Transition in a Strong Magnetic Background 
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We investigate tiie properties of tiie deconfining/ciiiral restoring transition for two flavor QCD in 
presence of a uniform background magnetic field. We adopt standard staggered fermions and a lattice 
spacing of the order of 0.3 fm. We explore different values of the bare quark mass, corresponding 
to pion masses in the range 200 - 480 MeV, and magnetic fields up to |e|i3 ~ 0.75 GeV . The 
deconfinement and chiral symmetry restoration temperatures remain compatible with each other 
and rise very slightly (< 2% for our largest magnetic field) as a function of the magnetic field. On 
the other hand, the transition seems to become sharper as the magnetic field increases. 

PACS numbers: 12.38.Aw, 11.15.Ha,12.38.Gc 



I. INTRODUCTION 

The study of strong interactions in presence of mag- 
netic background fields is relevant to many phenomeno- 
logical contexts. Large magnetic fields {B ^ 10^^ Tesla, 
i.e. \/]e\B ^ 1.5 GeV) may have been produced at 
the cosmological electroweak phase transition |l| and 
they may have influenced subsequent QCD transitions. 
Shghtly lower fields are expected to be produced in non- 
central heavy ion collisions, reaching up to 10"'^''^ Tesla at 
RHIC and up to - lO^^ Tesla at LHC l2,|3]. Large mag- 
netic fields, of the order of 10^° Tesla, are also expected 
in some neutron stars known as magnetars '4|. 

The influence of electric and magnetic fields on the chi- 
ral properties of the vacuum has been studied since some 
time, using various approximations or effective models of 
QCD |5l-ll2l|. predicting an enhancement of chiral symme- 
try breaking as a magnetic field is switched on. Recently 
new interesting phenomenology has been proposed, con- 
sisting in the appearance of an electric current parallel to 
the magnetic field in presence of deconfined quarks and 
local CP violations, induced e.g. by topological charge 
fluctuations [13l . [l4| . That is usually known as the chiral 
magnetic effect and experimental confirmations of it are 
currently being searched by heavy ion experiments [15| . 

An important issue is the infiuence of the magnetic 
field on the structure of the QCD phase diagram, in 
particular on the location and the nature of deconfine- 
ment and chiral symmetry restoration. Clarifying that 
in the case of strong magnetic fields is essential to cor- 
rectly predict the phenomenological consequences of the 
QCD transition on the evolution of the Universe dur- 
ing its early stages. Some computations exist, based on 
different approximations and QCD-like models [l6l - l22{ . 
which predict the possibility of a quite rich phenomenol- 
ogy, ranging from a possible splitting of deconfinement 
and chiral symmetry restoration to a sizable increase in 
the strength of the transition. However the various model 
predictions are not always consistent among themselves. 

A clarification of these issues may come from lat- 
tice QCD computations. A magnetic background field. 



contrary to an electric field or a finite baryon den- 
sity, does not give rise to technical difficulties such as 
a sign problem. The phase diagram in presence of a 
chromo-magnetic background field has been investigated 
in Rcfs. 23, 24]: the transition temperature decreases as 
a function of the external field, with deconfinement and 
chiral symmetry restoration staying strictly related to 
each other. Investigations in presence of (electro-) mag- 
netic fields have been done since long with the purpos e 
of studying the magnetic properties of hadrons [2^, ^^ , 
while some recent studies |27h30| have reported mostly 
on the chiral properties of the theory and about numeri- 
cal evidence for the chiral magnetic effect. 

In this paper we report on a first investigation of the 
QCD phase transition in presence of an (electro-) mag- 
netic background field. In order to do that, it is essen- 
tial to include dynamical quark contributions, since only 
quark fields, being electrically charged, are influenced by 
the magnetic field. We have considered Nf = 2 QCD 
with standard staggered fermions and different values of 
the quark masses, to appreciate how the effects of the 
magnetic field change as the mass spectrum changes (in 
the heavy quark limit the magnetic field becomes irrel- 
evant). In Section [n] we give some details about lattice 
QCD in presence of a background field and about our 
numerical setup. In Section Hill we present our numerical 
results and finally, in Section HVl we give our conclusions. 



II. NUMERICAL SETUP 

We consider Nf = 2 QCD, with quarks carrying differ- 
ent electric charges and coupled to a background (electro) 
magnetic (e.m.) field. The background field affects quark 
propagation and corresponds to a modification of the 
Dirac operator. In the continuum the covariant deriva- 
tive changes by inclusion of the e.m. A^^ field; on the lat- 
tice one has to add appropriate U{1) fields to the gauge 
link variables which parallel transport quarks fields from 
one lattice site to the other. In the case of a uniform 
magnetic field B, with different electric charges for the 



two flavors, g„ = 2[e[/3 and qd = — |e|/3 (|e| being the el- 
ementary charge) , the partition function of the (rooted) 
staggered fermion discretized version of the theory is 

Z(T, B)= f 23C/e"^'5 det Ali [B, q„] det M^ [B, qd] , (1) 

-. 4 

Mij[B, q] = amSij + - ^ rji^^ {u{B, q)i,M,^Stj-i, 

-u*{B,q\^,^,Ut,J,.,^o) ■ (2) 

TDIJ is the functional integration over the gauge link vari- 
ables C/„,^, Sg is the discretized pure gauge action (we 
consider a standard Wilson action) . The subscripts i and 
j refer to lattice sites, j> is a unit vector on the lattice 
and 7)1^1, are the staggered phases. Periodic (antiperi- 
odic) boundary conditions (b.c.) must be taken for gauge 
(fermion) fields along the Euclidean time direction, while 
spatial periodic b.c. are chosen for all fields. u{B, q)i^y 
are the gauge finks corresponding to the background U{\) 
magnetic field. We shall consider a constant magnetic 
field B ~ Bz and the following choice for the gauge field: 



.rl^ — tj'X 



Ap = for /i = .X, z, t 



(3) 



The corresponding C/(l) links on the lattice are: 



u{B,q)r 



^ia qBn^c 



u{B, q)n.^t = 1 for M = x, z, f (4) 



This choice corresponds to a magnetic flux a?B going 
through each plaquette in the x — y plane, except at the 
boundary (L^,, y, z, i), due to the periodic b.c. in the spa- 
tial directions. In order to guarantee the smoothness of 
the background field across the boundary and the gauge 
invariance of the fermion action the t/(l) gauge fields 
must be modified at the boundary of the x direction: 



"(i?,g)«,.=L. = 6-'"'*''^''"'' 



(5) 



and the magnetic field must be quantized, a?qB = 
2TTb/ LxLy, where b is an integer. That corresponds to 
taking the appropriate gauge invariant b.c. for fermion 
fields on the toru s I3l| (with the possible additional free 
phases 9x and 8y [31| set to zero). Given the two different 
values of g„ and qd, the quantization of B in our case is 
set by the d quark charge qd = — |e|/3. 



\e\B = BttT^ 



(6) 



T = l/{Nta) is the temperature and Lx = Ly = Lg. 

Our simulations have been carried out on 16^ x 4 lat- 
tices and for three different bare quark masses am = 
0.01335, 0.025 and 0.075. The corresponding (Goldstone) 
pion masses are arriT^ — 0.307(3), 0.417(3) and 0.707(3). 
The temperature T = 1/ [Nto) is changed by varying the 
lattice spacing through the inverse gauge coupling /3. 



Zero T estimates of the string tension, done at the same 
/3 values where the B — Q transition takes place, lead to 
a ranging from 0.29 to 0.31 fm as am is decreased, cor- 
responding to Tc{B — 0) ranging from 170 to 160 MeV. 
The corresponding values of the (Goldstone) pion mass 
are mjr k, 195, 275 and 480 MeV. For each quark mass we 
have done simulations using magnetic field correspond- 
ing to 6 = 0,8,16 and 24, i.e. for \e\B = 0,37rr2, 67rT2 
and QttT^ respectively. Thus, for the lightest pion mass 
our magnetic field reaches values up to \e\B ss 19 m^, i.e. 
■y/]e|S « 850 MeV in physical units. Note that, since we 
are working with a fixed value oi Nt, B changes while 
changing temperature as a~^ oc T^. However, for all the 
quark masses the range of couplings (hence of a) that we 
explore corresponds to a < 2% change in T and hence the 
magnetic field only changes at most by a few percent. 

The Rational Hybrid Monte-Carlo algorithm has been 
used to simulate rooted staggered fermions: we need to 
treat separately each fiavor and thus take the fourth root 
of the fermion determinant. Typical statistics are of the 
order of 10k molecular dynamics trajectories. 



III. NUMERICAL RESULTS 

In Figs. midland [3] we show the behavior of (■(/'V') (aver- 
age of the u and d quark condensates) and of the Polyakov 
loop for different magnetic fields and am. — 0.075, 0.025 
and 0.01335. Results are presented as a function of the 
inverse gauge coupling /3: we recall that T is an increas- 
ing function of /3, a conversion into physical units will be 
presented later. For all temperatures ('0V') increases as 
a function oi B, as expected from analytic predictions. 
Interpreting the drop of the condensate as the signal for 
chiral symmetry restoration, we infer that the transition 
temperature increases as a function of B. A sharper drop 
is observed at the highest fields explored, especially for 
the lowest quark masses, indicating a sizable increase in 
the strength of the transition; that is visible from the be- 
havior of the disconnected chiral susceptibility (Fig. 01). 

The Polyakov loop P is a pure gauge quantity, coupled 
to the magnetic field only through quark loops, hence its 
behavior is less trivial to predict. From Figs. [1] [2] and [3] 
we see that, while at low T it decreases as a function of 
B (as one would expect qualitatively from the fact that 
(■00) increases), at high T it increases. Such behavior is 
qualitatively similar to what obtained in Ref. [20| by a 
PNJL model analysis and should be further investigated, 
e.g. by determining the renormalized Polyakov loop. If 
we interpret the rise of P as the onset of deconfinement, 
we infer that the shift and the increase in strength of 
the deconfinement transition is similar to what observed 
for the chiral transition. Data obtained for the Polyakov 
loop susceptibility lead to similar conclusions (see Fig. [5]). 

In Table H] we report the pseudocritical couplings /3c for 
deconfinement and chiral restoration obtained by fitting 
the peak of the susceptibilities by a quadratic function. 
We have also determined /3c looking for the inflection 
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FIG. 1: Chiral condensate and Polyakov loop for am = 0.075 pjg 4. Disconnected (^V) susceptibility for am = 0.01335 
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FIG. 2: Same as in Fig. [T] for am = 0.025 
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FIG. 3: Same as in Fig.[T]for am = 0.01335 



point of observables, by means of polynomial fits, obtain- 
ing compatible results within errors. Data obtained for /3c 



FIG. 5: Polyakov loop susceptibility for am, = 0.01335 



confirm that no appreciable separation of chiral restora- 
tion and deconfinement is induced by the the background 
field, at least for the explored field strengths. 

From the values of /3c we obtain the ratio Tc{B) /Tc{0) 
as a function of the dimensionless ratio eB/T^, as re- 
ported in Fig. [5] the 2-loop /3-function has been used for 
the conversion. A direct determination of the physical 
scale on T = lattices is preferable but would require 
very precise measurements to appreciate T variations of 
the order of percent. On the other hand, given the small 
scale range explored, the approximation is reasonable and 
no qualitative change is expected. 

Fig. |6] shows that the change in Tc is small and of 
the order of a few percent at the highest explored fields. 
Moreover, there seems to be a saturation as the chi- 
ral limit is approached: results for am — 0.01335 and 
am = 0.025 stay onto each other. Notice that this is true 
if we plot results as a function of \e\B/T'^: had we used 
\e\B /m^ results at different masses would have been dif- 
ferent: the highest B is about 20 m^ for am, = 0.01335 
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TABLE I: Pseudocritical couplings obtained by fitting the 
peak of the chiral condensate or Polyakov loop susceptibilities. 
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FIG. 6; Tc{B) for different quark masses. The solid curve is 
a power law fit to the lightest quark data (see text). 




FIG. 7: Reweighted plaquette distribution at /3c as a function 
of the external field at am — 0.01335 on a 16"^ x 4 lattice. 



and 10 mj for am = 0.025. This suggests that, at least 
for the strong fields and for the pion masses explored, the 
relevant scale governing the effect of the magnetic field 
on the shift of the transition is T itself and not ttLt^ . 
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FIG. 8: Monte-Carlo histories of the plaquette at 3 different 
/3 values around the transition for 6 = 24 and am = 0.01335. 



Trying to understand the dependence of Tc on B, we 
have fitted our data for am = 0.01335 according to 



TciO) 



l + A 



\e\B 

rp2 



(7) 



finding that a = 1.45(20) and A - 1.3 10"^. 

Finally we discuss about the nature of the transition. 
At B = it is still unclear if a weak first order transition 
is present in the chiral limit [32|, |33|, however no clear 
signal for a finite latent heat has been found on avail- 
able lattice sizes, hence the first order transition, even if 
present, is so weak to be of poor phcnomcnological rele- 
vance. On the other hand our results show that the intro- 
duction of a magnetic field makes the transition sharper. 
The question is if large fields can turn the transition into 
a first order strong enough to be clearly detectable. 

To that aim we have analyzed the reweighted plaque- 
tte distribution at the critical couplings and for different 
values of B: results are shown in Fig. [5] The single peak 
distribution, which is present at zero or small magnetic 
field, turns into a double peak distribution, typical of a 
first order transition, for the largest B explored; also the 
Monte-Carlo histories of the plaquette. Fig. [SJ present 
signals of a metastable behavior. We can consider that 
as an indication but not as a final answer: numerical sim- 
ulations on larger lattice sizes are necessary to clarify if 
the double peak structure survives the thermodynamical 
limit and for a proper finite size scaling analysis. 



IV. CONCLUSIONS 

We have presented results from an investigation of 
the Nf = 2 QCD phase diagram in presence of a mag- 
netic background field. We have explored different quark 
masses, corresponding to m,r ranging from 200 MeV to 
480 MeV, and different magnetic fields, with y/\e\B up 
to about 850 MeV (|e|B - 20 to^ for the lightest mass). 



Main results can be summarized as follows: the transi- 
tion temperature increases slightly ( < 2% at the highest 
field) and no evidence is found, within the range of ex- 
plored fields, for a disentanglement of chiral symmetry 
restoration and deconfinement. Tc{B)/Tc{0) as a func- 
tion of \e\B/T'^ shows negligible dependence on TOtt for 
the two lowest masses, and is well described by a power 
law Tc{B)/TciO) = 1 + A{\e\B/T^)°' with a - 1.45(20). 
The transition becomes sharper with some preliminary 
evidence for a first order transition, in the form of double 
peak distributions, at the highest fields explored: such in- 
dications should be clarified by future studies on larger 
spatial volumes and by a finite size scaling analysis. 

Regarding the comparison with model predictions, our 
results show partial agreement with some of the results 
reported in Ref. [21| and in Ref . 20] : the deconfinement 
and chiral restoring temperatures both increase, even if 
we do not see any sign for a faster grow and splitting of 
the chiral transition till \e\B ~ 20 m^. Also the observed 
increase in the strength of the transition is common to 
some models [iE;[2l|. We stress the qualitatively different 
behavior which is observed in numerical simulations with 



a background chromo- magnetic field, where Tc decreases 
as a function of the external field [2j, [2J] . 

Our results have been obtained using a standard stag- 
gered discretization and a coarse lattice, with a lattice 
spacing a ^ 0.3 fm. Apart from possible systematic ef- 
fects related to the fourth root trick, fiavour breaking 
discretization effects may play an important role, with 
a distorted hadron spectrum that could partially modify 
the effect of the magnetic field on QCD thermodynamics. 
For this reason it will be important to confirm our results 
in the future by using different lattice discretizations. 
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